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Abstract 

This article concerns a phenomenon of elementary quantum mechanics that 
is quite counter-intuitive, very non-classical, and apparently not widely known: 
a quantum particle can get reflected at a potential step downwards. In contrast, 
classical particles get reflected only at upward steps. As a consequence, a quantum 
particle can be trapped for a long time (though not forever) in a region surrounded 
by downward potential steps, that is, on a plateau. Said succinctly, a quantum 
particle tends not to fall off a table. The conditions for this effect are that the wave 
length is much greater than the width of the potential step and the kinetic energy 
of the particle is much smaller than the depth of the potential step. We point out 
how the topic is accessible with elementary methods, but also with mathematical 
rigor and numerically. 

PACS: 03.65.-w, 03.65.Nk, 01.30.Rr. Key words: Schrodinger equation; potential 
step; confining potential; reflection and transmission coefficients. 

1 Introduction 

Suppose a quantum particle moves towards a sudden drop of potential as in Figure [[J 
with the particle arriving from the left. As a classical analogue, we may imagine a ball 
rolling towards the edge of a table. Will it fall down or be reflected? Classically, the 
ball is certain to fall down, but a quantum particle has a chance to be reflected. That 
sounds paradoxical because the particle turns around and returns to the left under a 
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force pointing to the right! Under suitable conditions reflection even becomes close to 
certain. This non-classical, counter-intuitive quantum phenomenon we call "paradoxical 
reflection," or, when a region is surrounded by downward potential steps, "paradoxical 
confinement" — where "paradoxical" is understood in the sense of "counter-intuitive," 
not "illogical." It can be derived easily using the following simple reasoning. 
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Figure 1: A potential V(x) containing a step downwards 

Suppose the particle moves in 1 dimension, and the potential is a rectangular step 
as in Figure [U 

V(x) = -AEG(x) (1) 

with O the Heaviside function and AE > 0. A wave packet coming from the left 
gets partially reflected at the step and partially transmitted. The size of the reflected 
and the transmitted packets can be determined by a standard textbook method of 
calculation (e.g., [HE]), the stationary analysis, replacing the wave packet by a plane 
wave of energy E and solving the stationary Schrodinger equation. The transmitted 
and reflected probability currents, divided by the incoming current, yield the reflection 
and transmission coefficients R > and T > with R + T = 1. We give the results 
in Section [2] and observe two things: First, R ^ 0, implying that partial reflection 
occurs although the potential step is downwards. Second, R even converges to 1, so that 
reflection becomes nearly certain, as the ratio Ej AE goes to zero. Thus, paradoxical 
reflection can be made arbitrarily strong by a suitable choice of parameters (e.g., for 
sufficiently big AE if E is kept fixed). 

If it sounds incredible that a particle can be repelled by a potential step downwards, 
the following fact may add to the amazement. As derived in [2j p. 76], the reflection 
coefficient does not depend on whether the incoming wave comes from the left or from 
the right (provided the total energy and the potential are not changed). Thus, a step 
downwards yields the same reflection coefficient as a step upwards. (But keep in mind 
the difference between a step upwards and a step downwards that at a step upwards, 
also energies below the height of the step are possible for the incoming particle, a case 
in which reflection is certain, R = 1.) 

To provide some perspective, it may be worthwhile to point to some parallels with 
quantum tunneling: there, the probability of a quantum particle passing through a 
potential barrier is positive even in cases in which this is impossible for a classical 
particle. In fact, paradoxical reflection is somewhat similar to what could be called 
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anti-tunneling, the effect that a quantum particle can have positive probability of being 
reflected by a barrier so small that a classical particle would be certain to cross it. 

In the rest of this article, we address the following questions: Is paradoxical reflection 
a real physical phenomenon or an artifact of mathematical over-simplification? (We 
will look at numerical and rigorous mathematical results.) How does it depend on 
the parameters of the situation: the width L (see Figure [2]) and the depth AE of 
the potential step, the wave length A and the width a of the incoming wave packet? 
Why does this phenomenon not occur in the classical regime? That is, how can it be 
that classical mechanics is a limit of quantum mechanics if paradoxical reflection occurs 
in the latter but not the former? And, could one use this phenomenon in principle 
for constructing a particle trap? In spring 2005, these questions gave rise to lively and 
contentious discussions between a number of physics researchers visiting the Institut des 
Hautes Etudes Scientifiques near Paris, France; these discussions inspired the present 
article. 



2 Stationary Analysis of the Rectangular Step 

We begin by providing more detail about the stationary analysis of the rectangular step 
(PQ), considering the time- independent Schrodinger equation (m = mass) 

E1>{x) = -Jtl>"(x) + V(x)1>(x). (2) 

This can be solved in a standard way: for x < 0, let if) be a superposition of an incoming 
wave e lklX and a reflected wave Be~ lklX , while for x > 0, let ip be a transmitted wave 
Ae lk2X , with a possibly different wave number k 2 . Indeed, from ([2]) we obtain that 

k ± = V2mE/h, k 2 = \/2m(E + AE)/h. (3) 

The value E > is the kinetic energy associated with the incoming wave. The coef- 
ficients A and B are determined by continuity of ip and its derivative if)' at x = to 
be 

ki + k 2 ki + k 2 

The reflection and transmission coefficients R and T are defined as the quotient of the 
quantum probability current j = {h/m)\m{%p*il)') associated with the reflected respec- 
tively transmitted wave divided by the current associated with the incoming wave, 

R= \jlA^ T = 3 -^. (5) 

Jin Jin 

Noting that j t ra = hk 2 \A\ 2 /m, j refl = —hki\B\ 2 /m, j in = hki/m, we find that 

R= \B\ 2 = 1-^\A\\ T = ^\A\ 2 . (6) 
ki ki 
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Note that both R and T lie in the interval [0, 1], and that R + T — 1. By inserting (T4J) 
into (jSJ), we obtain 

(h + k 2 f - Ahk 2 = (fc 2 - fcQ 2 

(A* + fc 2 ) 2 (fex + k 2 f 1 j 

and make two observations: First, R ^ 0, implying that reflection occurs, if k\ ^ k 2 , 
which is the case as soon as AE ^ 0. Second, R even converges to 1, so that reflection 
becomes nearly certain, as the ratio r := Ej AE tends to zero; that is because 



R 
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k 2 + h I \ y/E + AE + VE \ Vr + l + ^fr 



since both the numerator and the denominator tend to 1 as r — > 0. This is the simplest 
derivation of paradoxical reflection. 



3 Soft Step 

For a deeper analysis of the effect, we will gradually consider increasingly more realistic 
models. In this section, we consider a soft (or smooth, i.e., differentiable) potential step, 
as in Figure [21 for which the drop in the potential is not infinitely rapid but takes place 
over some distance L. The result will be that paradoxical reflection exists also for soft 
steps, so that the effect is not just a curious feature of rectangular steps (which could not 
be expected to ever occur in nature). Another result concerns how the effect depends 
on the width L of the step. 
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Figure 2: A potential containing a soft step 

To study this case it is useful to consider the explicit function 

AE / x \ 

V(x) = -— (l+tanh-J, (9) 

depicted in Figure [3 (Recall that tanh = sinh / cosh converges to ±1 as x — > ±oo.) The 
reflection coefficient for this potential can be calculated again by a stationary analysis, 
obtaining from the time-independent Schrodinger equation (j2J) solutions ip(x) which are 
asymptotic to e lfclX + Be~ lklX as x — ► — oo and asymptotic to Ae lk2X as x — > oo, i.e., 
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lim a ._ >oo (-0(a;) — Ae lk2X ) = 0. The calculation is done in [21 p. 78]: The values k± and k 2 
are again given by ([3]), and the reflection coefficient turns out to be 

^sinh(f(A; 2 + A;i)L) J 1 ' 

From this and ([3]) we can read off that, again, R ^ for AE ^ 0, and i? — > 1 as S — > 
while AS and L are fixed (since then /c! — > 0, fc 2 — > v2rnAE /h, so both the numerator 
and the denominator tend to sinh( | a/ ImAEL /h)). As AS — > oo while S and L are 
fixed, -> exp(-2vr V2mEL/h) because for large arguments sinh m ^ exp. 

In addition, we can keep E and AE fixed and see how R varies with L: In the limit 
L — > 0, ffTUl) converges to ([7]) because sinh(aL) w aL for L 1 and fixed a; this is 
what one would expect when the step becomes sharper and ([HD converges to ([T]). In the 
limit L — > oo, R converges to because for fixed (3 > a > 

sinh(aL) e aL - e~ aL e ^~^ L - ^~ a ~P) L 

= = > fill 

sinh(/?L) eP L -e-0 L 1 - e~^ L ' 1 ' 

as the numerator tends to and the denominator to 1. Thus, paradoxical reflection 
disappears for large L; in other words, it is crucial for the effect that the drop in the 
potential is sudden. 

Moreover, ( TlUj) is a decreasing function of L, which means that reflection will be the 
more probable the more sudden the drop in the potential is. To see this, let us check 
that for (3 > a > and L > the function f(L) = sinh(aL)/ sinh(/3L) is decreasing: 

df a cosh(aL) sinh(/?L) — (3 sinh(aL) cosh(/3L) 



dL smh z ((3L) 



because 



° <r4^. (13) 



tanh a tanh (3 

as xj tanhx is increasing for x > 0. 

How about soft steps with other shapes than that of the tanh function? Suppose 
that the potential V(x) is a continuous, monotonically decreasing function such that 
V(x) — > as x — > — oo and V(x) — > —AE as x — > +oo. To begin with, we note 
that the fact, mentioned in the introduction, that the reflection coefficient is the same 
for particles coming from the left or from the right, still holds true for such a general 
potential [21 p. 76]. This suggests that paradoxical reflection occurs also for general 
potential steps. Unfortunately, we do not know of any general result on lower bounds 
for the reflection coefficient R that could be used to establish paradoxical reflection in 
this generality. However, an upper bound is known [51 eq. (82)], according to which R 
is less than or equal to the reflection coefficient d7j) of the rectangular step. This agrees 
with our observation in the previous paragraph that reflection is the more likely the 
sharper the step. 
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4 Wave Packets 



Another respect in which we can be more realistic is by admitting that the wave function 
with which a quantum particle reaches a potential step is not an infinitely-extended 
plane wave e lklX but in fact a wave packet of finite width er, for example a Gaussian 
wave packet 

^■ m {x) = G^{x) l l 2 e^ x (14) 
with Gfj„cr the Gauss function with mean \x and variance a 2 , 

G^{x) = -^ e -(z-rf/^ 2 . (15) 

V27TO" 



Suppose this packet arrives from the left and evolves in the potential V(x) according to 
the time-dependent Schrodinger equation 

ih^(x, t) = ~P^{x, t) + V(x) rf,(x, t) . (16) 
at 2m ox z 

Ultimately, as t — > oo, there will be a reflected packet iprefi in the region x < moving 
to the left and a transmitted packet iptr& in the region x > moving to the right, and 
thus the reflection and transmission probabilities are 

J R=||^cfi|| 2 , T=||^ tra || 2 , (17) 

with |H| 2 = n>(:r)| 2 d*. 

Because of the paradoxical feel of paradoxical reflection, one might suspect at first 
that the effect does not exist for wave packets but is merely an artifact of the stationary 
analysis. We thus address, in this section, the question as to how wave packets behave, 
and whether the reflection probability (JTTj) agrees with the reflection coefficient discussed 
earlier. We begin with the numerical evidence confirming paradoxical reflection. 



4.1 Numerical Simulation 

A numerical simulation of a wave packet partly reflected from a (soft) step downwards 
is shown in Figure [3j The simulation starts with a Gaussian wave packet moving to 
the right and initially located on the left of the potential step. After passing the step, 
there remain two wave packets, no longer of exactly Gaussian shape, one continuing to 
move to the right and the other, reflected one returning to the left. For the choice of 
parameters in this simulation, the transmitted and reflected packet are of comparable 
size, thus providing evidence that there can be a substantial probability of reflection 
at a step downwards (even for wave packets of finite width). That is, the numerical 
simulation confirms the prediction of the stationary analysis. 
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Figure 3: Numerical simulation of the time- dependent Schrodinger equation for the soft 
step potential (jUJ). The picture shows ten snapshots of \ip\ 2 (black lines) at different 
times before, during, and after passing the potential step. (Order: left column top 
to bottom, then right column top to bottom.) It can clearly be seen that there is a 
transmitted wave packet and, "paradoxically," a reflected wave packet. The initial wave 
function is a Gaussian wave packet centered at x = 0.1 with o = 0.01 and ko = 2007T. 
The simulation assumes infinite walls at x = and x = 1. The parameters are m = 1, 
L = 1/100 and AE = 18.4 E, and the linear mesh has iV = 1000 points. The additional 
lines in the figures depict the potential in arbitrary units. The figure shows the modulus 
of the wave packet at times: 1, 4, 8, 11, 12, 13, 16, 18, 22, 27 At in arbitrary time units. 

4.2 But Is It for Real? 

We now point out how rigorous mathematics confirms paradoxical reflection as a conse- 
quence of the Schrodinger equation. We thus exclude the possibility that it was merely 
numerical error that led to the appearance of paradoxical reflection for wave packets. 

Do not think the worry that numerical errors may lead to the wrong behavior of a 
wave packet was paranoid: There are cases in which exactly this happens. Here is an 
example involving the potential V(x) = —x 2 : In some numerical simulations, a wave 
packet starting near the origin and moving to the right gets reflected at some point and 
returns to the origin, whereas in the rigorous solution of the Schrodinger equation the 
wave packet never returns, but accelerates more and more. As illustrated in Figure HI the 
numerically predicted time of reflection depends on the mesh width (i.e., the numerical 
resolution): the higher the resolution, the later the predicted reflection. 
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Figure 4: An example of how numerical error may lead to wrong predictions. The plot 
shows numerical results about the evolution of the average position of a Gaussian wave 
packet (initially as in Figure EJ) with the parabolic potential V(x) = —50k^(x — 0.3) 2 
as a function of time for several numbers of nodes of the linear mesh, N. The black 
dashed line corresponds to the evolution for V(x) = 0. We see that there is a spurious 
reflection of the wave packet due to the finite size effects of the numerical simulation. 
Note that in the figure only in the iV = 10, 000 case the wave function collides with the 
x = 1 wall. In the other cases the reflection occurs far from the x = 1 wall (see for 
instance the N = 500 case). 

We return to the mathematics of paradoxical reflection. The rigorous mathematical 
analysis of scattering problems of this type is a fairly complex and subtle topic. The 
main techniques and results (also for higher dimensional problems) are described in [HE], 
and the mathematical results relevant to potentials of the step type can be found in [6]. 
The reflection probability R of eq. (Tl~7l) is given in terms of the plane wave reflection 
coefficients R(k\) by the following formula, expressing exactly what one would intuitively 
expect: 

poo 

R= / dfcii^OlJtofa)! 2 . (18) 
Jo 

The same formula holds with all R's replaced by T's. In (fT8l) . R(ki) is given by the 
stationary analysis, as in ([7j) or fflUj) . with k% expressed in terms of k\ and AE, = 
+ 2mAE/h 2 ; and ^- m (ki) is the Fourier transform of the incoming wave packet 

^ in (x). 
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(In brackets: To be precise, the incoming packet ip in (x,t) is defined as the free 
asymptote of ip(x,t) for t — > — oo, i.e., ip in (x,t) evolves without the potential, 

ih ^L( x t ) = - — ^(x t) (19) 
n dt 1 ' J 2m Ox 2 [ ' J ' 1 J 

and 

lim \\^ n {;t)-^{;t)\\ =0. (20) 

t— >— oo 

Similarly, ^ re fl + is the free asymptote of ip f° r t — > +oo. When we write ^ in (x), 
we mean to set t = 0; note, however, that ffT8l) actually does not depend on the t- value 
as, by the free Schrodinger equation (fl9l) . ^i n (fc, t) = exp(— ithk 2 /2m) ipi n (k, 0) and thus 
|V>in(fc, ^)P = iV'in(fc) 0)| 2 . Since we assumed that the incoming wave packet comes from 
the left, ipi n is a "right-moving" wave packet consisting only of Fourier components with 
k > 0.) 

From (j!8p we can read off the following: If the incoming wave packet consists only 
of Fourier components k\ for which R(k\) > 1 — e for some (small) e > 0, then also 
R > 1 — e. More generally, if the incoming wave packet consists mainly of Fourier 
components with R(ki) > 1 — e, that is, if the proportion of Fourier components with 
R(h) > 1 - e is 

poo 

/ dAi9(%) - (1 - e)) $Uh)\ 2 = 1-5, (21) 
Jo 

then R > 1 — e — 5 because 

poo poo 

/ dfcii2(fci)fc(*i)| 2 > / dk 1 R(k 1 )e(R(k 1 )-(l-s))\^ m (k 1 )\ 2 > 
Jo Jo 

/■oo 

> / dh (1 - e) e(R(h) - (1 - 5)) |4(^i)| 2 = (1 - - 5) > 1 - £ - S . 
Jo 

Therefore, whenever the stationary analysis predicts paradoxical reflection for cer- 
tain parameters and values of k\, then also wave packets consisting of such Fourier 
components will be subject to paradoxical reflection. 



5 Parameter Dependence 

Let us summarize and be explicit about how the reflection probability R from a potential 
step downwards depends on the parameters of the situation: the mean wave number k\ 
and the width a of the incoming wave packet, and the depth AE and width L of the 
potential step. We claim that R is close to 1 in the parameter region with 

— > L (22a) 
k\ 

h 2 k 2 

AE^>^ = E (22b) 
2m 

a > . (22c) 
h 
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Note that l/k\ is (up to the factor 2tt) the (mean) wave length A. 
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Figure 5: The region (shaded) in the plane of the parameters u and v, defined in (|23|) . 
in which the reflection probability (124"|) exceeds 99 percent. The horizontally shaded 
subset is the region in which condition (1251) holds. 

To derive this claim from ffTDT) and ffTS]) , consider first the case a — > oo of a very very 
wide packet. For such a packet, its Fourier transform is very sharply peaked at k±. The 
reflection coefficient R given by (flOl) depends on the parameters k\,L, AE,m only in 
the dimensionless combinations 



u 



v = %y/2mAEL/h. 



that is, 



R = R(u,v) 



sinh(v / w 2 + v 2 — u) 



(23) 



(24) 



sinh(v / ^ 2 + v 2 + u) 

Figure [5] shows the region in the uv plane in which R > 0.99. As one can read off from 
the figure, for (u,v) to lie in that region, it is sufficient, for example, that 

u < 10~ 3 and v > 10 3 u . (25) 

More generally, for R(u, v) to be very close to 1 it is sufficient that m < 1 and v ^> u, 
which means (I22al) and f!22bj) . To see this, note that 



sinh(v / u 2 + v 2 — u) 



sinh("\/w 2 + v 2 + u — 2u) 



so that 



U, V 



smh(Vu 2 + v 2 + u) cosh(2u) — cosh(-\/ u 2 + v 2 + u) sinh(2w) 
sinh(v / w 2 + v 2 — u) , sinh(2w) 



(26) 



smh(\/u 2 + v 2 + u) 



cosh (2m) — 



tanh(\/M 2 + v 2 + u) 



(27) 
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Suppose that m < 1. Then Taylor expansion to first order in u yields 



y/R{u,v) « 1 - -\- . (28) 
tanhf 

If v is of order 1, this is close to 1 because u < 1. If, however, v is small, then tanhw 
is of order v, and the right hand side of (|28p is close to 1 when u/v <C 1. Thus, when 
(I22al) and (122bl) are satisfied \/R is close to 1 , and thus so is R. 

Now consider a wave packet that is less sharply peaked in the momentum representa- 
tion. If it has width a in position space then, by the Heisenberg uncertainty relation, it 
has width of order 1/a in Fourier space. For the reflection probability to be close to one, 
the wave packet should consist almost exclusively of Fourier modes that have reflection 
coefficient close to one. Thus, every wave number k\ in the interval, say, [k\ — - , h\-\- —} 
should satisfy (122aj) and (I22bl) . This will be the case if — is small compared to ki, or 
a ^> 1/ki. Thus, ( I22cp . which is merely what is required for (|14p to be a good wave 
packet, i.e., an approximate plane wave, is a natural condition on a for keeping R close 
to 1. 



6 The Classical Limit 

If paradoxical reflection exists, then why do we not see it in the classical limit? On 
the basis of (I22p we can understand why: Classical mechanics is a good approximation 
to quantum mechanics in the regime in which a wave packet moves in a potential that 
varies very slowly in space, so that the force varies appreciably only over distances much 
larger than the wave length. For paradoxical reflection, in contrast, it is essential that 
the length scale of the drop in the potential be smaller than the wave length. For further 
discussion of the classical limit of quantum mechanics, see [7] . 



7 A Plateau as a Trap 

Given that a quantum particle will likely be reflected from a suitable potential step 
downwards, it is obvious that it could be trapped, more or less, in a region surrounded 
by such a potential step. In other words, also potential plateaus, not only potential 
valleys, can be confining. To explore this possibility of "paradoxical confinement," we 
now consider a potential plateau 

V(x) = -AE(e(x-a) + e(-x-a)), (29) 

depicted in Figure El 

A particle starting on the plateau could remain there — at least with high probability — 
for a very long time, much longer than the maximal time r c i that a classical particle 
with energy E would remain on the plateau, which is 

[2m , . 

r d = ay ~y , (30) 
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Figure 6: Potential plateau 



independently of the height of the plateau. Instead, the time the quantum particle likely 
remains on the plateau is of the order a/ AE/E t c \ and is thus much larger than t c \ if 
the height AE is large enough. One can say that a quantum particle on a table tends 
not to fall down. In this formulation, the "table" is supposed to be of an appreciable 
height, i.e., AE ^> E. In fact, as we shall argue in the subsequent sections, a quantum 
particle starting on the plateau should leave the plateau at the rate r q ^ with the decay 
time 

V2mA£ 1 [AE 



T, 



qu 



T d . (31) 

AE 4 V E cl K J 

The lifetime (13T1) can be obtained in the following semi-classical way: Imagine a 
particle traveling along the plateau with the speed ^/2E/m classically corresponding to 
energy E, getting reflected at the edge with probability R given by ([7]), traveling back 
with the same speed, getting reflected at the other edge with probability R, and so on. 
Since the transmission probability T = 1 — R corresponding to ([7]) is 



A^jE/AE + higher powers of Ej AE . 



(32) 



a number of reflections of order {^J E / AE)^ 1 should typically be required before trans- 
mission occurs, in qualitative agreement with fl3TT) . In fact, the transmission probability 
of T = 4a/ E/AE, when small, corresponds to a decay rate T/r c \ and hence to the decay 
time t c \/T given by ( 13TI) . 



V(x)\ 



AE 



Figure 7: Potential well 

One must be careful with this reasoning, since applied carelessly it would lead to the 
same lifetime for the potential well depicted in Figure [7J as for the potential plateau. 
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(That is because the reflection probability at a potential step upwards is, as already 
mentioned, the same as that at a potential step downwards.) However, the potential well 
possesses bound states for which the lifetime is infinite. In this regard it is important 
to bear in mind that the symmetry in the reflection coefficient derived in [2] always 
involves incoming waves at the same total energy E > 0; for a potential well it would 
thus say nothing about bound states, which have E < 0. 

This is a basic difference between confinement in a potential well and paradoxical 
confinement on a potential plateau: In the well, the particle has positive probability 
to stay forever. Mathematically speaking, the potential well has bound states (i.e., 
eigenfunctions in the Hilbert space L 2 (M) of square-integrable functions), whereas the 
potential plateau does not. For the potential well, the initial wave packet will typically 
be a superposition ip = Abound + ^scattering of a bound state (a superposition of one or 
more square-integrable eigenfunctions) and a scattering state (orthogonal to all bound 
states); then || Abound || 2 is the probability that the particle remains in (a neighborhood of) 
the well forever. In contrast, because of paradoxical reflection, the potential plateau has 
metastable states, which remain on the plateau for a long time but not forever. Indeed, 
we will show that (1311) is the lifetime. Our tool for the proof will be a method similar to 
the stationary analysis of Section 2, using special states lying outside the Hilbert space 
L 2 (M) (as do the stationary states of Section 2). And again like the stationary states of 
Section 2, the special states are similar to eigenfunctions of the Hamiltonian: they are 
solutions of the time-independent Schrodinger equation (j2]), but with complex "energy"! 

8 Eigenfunctions with "Complex Energy" 

We now derive the formula fl3"TT) for the lifetime r = r qu from the behavior of solutions 
to the eigenvalue equation (j2J), but with complex eigenvalues. To avoid confusion, let us 
now call the eigenvalue Z instead of E; thus, the equation reads 

Z^(x) = -g^"(x) + V(x) 1>(x) , (33) 

where V is the plateau potential as in (f29|) . Such "eigenfunctions of complex energy" 
were first considered by Gamow [HI E] for the theoretical treatment of radioactive alpha 
decay. 

The fact that the eigenvalue is complex may be confusing at first, since the Hamilto- 
nian is a self-adjoint operator, and it is a known fact that the eigenvalues of a self-adjoint 
operator are real. However, in the standard mathematical terminology for self-adjoint 
operators in Hilbert spaces, the words "eigenvalue" and "eigenfunction" are reserved for 
such solutions of ( |33i) that if> is square-integrable (= normalizable), i.e., ip £ L 2 (M). In 
this sense, all eigenvalues must be real indeed; for us this means that any solution if) 
of (|33|) for Z G C \ K is not square-integrable. In fact, even the eigenfunctions with 
real eigenvalue E considered in (j2j) were not square-integrable, which means that they 
do not count as "eigenfunctions" in the mathematical terminology, and do not make 
the number E an "eigenvalue." Instead, E is called an element of the spectrum of the 
Hamiltonian. Still, the spectrum of any self-adjoint operator consists of real numbers, 
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and thus ZgC\1 cannot belong to the spectrum of the Hamiltonian. So, the eigen- 
values Z we are talking about are neither eigenvalues in the standard sense, nor even 
elements of the spectrum. (Nevertheless we continue calling them "eigenvalues," as they 
satisfy (|33|) for some nonzero function.) 

Let us explain how these complex eigenvalues can be useful in describing the time 
evolution of wave functions. Consider an eigenfunction ip with a complex eigenvalue Z. 
It generates a solution to the time-dependent Schrodinger equation by defining 

tjj(x,t) =e- iZt/h ^(x,0). (34) 

The function grows or shrinks exponentially with time, with rate given by the imaginary 
part of Z. More precisely, 

\^(x,t)\ 2 = e 21raZt/h \ij(x,0)\ 2 } (35) 

so that 2lmZ/h is the rate of growth of the density \ip(x,t)\ 2 . For those eigenf unctions 
relevant to our purposes, the imaginary part of Z is always negative, so that ip shrinks 
with time. In particular, the amount of \ijj\ 2 on top of the plateau decays with the expo- 
nential factor that occurs in fl35l) . Assuming that \ip\ 2 is proportional to the probability 
density at least in some region around the plateau (though not on the entire real line) 
for a sufficiently long time, and using that the lifetime r for which the particle remains 
on the plateau is reciprocal to the decay rate of the amount of probability on top of the 
plateau, we have that 

T = -2i£z- < 36 > 

From (134p we can further read off that the phase of ip(x,t) at any fixed x rotates 
with frequency KeZ/h, while for eigenf unctions with real eigenvalue E it does so with 
frequency E/h, which motivates us to call ReZ the energy and denote it by E. Thus, 

Z — E — i— . (37) 

Below, we will determine r by determining the relevant eigenvalues Z, i.e., those corre- 
sponding to decay eigenfunctions, see below (HUj) . 

How can it be that \ip\ 2 shrinks everywhere? Does that not conflict with the conser- 
vation of the total probability? Apart from the fact that \ip\ 2 cannot be the probability 
density anyway since ip is not square-integrable, the shrinking of ip means that there 
is a flow of \ip | 2 to infinity (that is not compensated by an equally large flow of \i/j\ 2 
from infinity). This is indeed similar to the situation we want to consider, in which the 
amount of \ijj\ 2 on top of the plateau continuously shrinks due to a flow of \ip\ 2 away from 
the plateau. It also suggests that ip should grow exponentially in space as x — > ±oo: 
The density at great distance from the plateau would be expected to agree with the flow 
off the plateau in the distant past, which was exponentially larger than in the present 
if the wave function on top of the plateau shrinks exponentially with time. As we will 
see now, the eigenfunctions behave exactly that way. 



14 



The general solution of fl33|) with Z E C except Z = or Z = —AE ia3 

{B_e~ lkx + C_e lkx when x < —a, 
A + e lkx + A_e~ lkx when — a < x < a, (38) 
B + e lkx + C + e~ lkx when x > a, 

where 

k = V2mZ/h and k = \/2m(Z + AE)/h (39) 

with the following (usual) definition of the complex square root: Given a complex num- 
ber ( other than one that is real and < 0, let a/C denote the square root with positive 
real part, Re\/( > 0. For ( < 0, we let = iyjCl- (Since (138]) remains invariant 
under changes in the signs of k and k, choosing the positive branch for the square roots 
is not a restriction for the solutions.) We are interested only in those solutions for Z 
with ReZJ = E > 0; these are the ones that should be relevant to the behavior of states 
starting out on the plateau (with positive energy). Nevertheless, for mathematical sim- 
plicity, we will also allow Re Z < but exclude any Z that is real and negative or zero. 
For any Z G C \ (— oo,0] we have that Kek > 0, so that the probability current j 
associated with exp(ikx) is positive (i.e., a vector pointing to the right). Since we do 
not want to consider any contribution with a current from infinity to the plateau, we 
assume that 

C + = C- = . (40) 

We thus define a decay eigenfunction or Gamow eigenfunction to be a nonzero function 
ip of the form (I3"8"|) with (1351 and C± = 0, satisfying the eigenvalue equation (1331) except 
at x = ±a (where ip" does not exist) for some Z G C \ (— oo,0], such that both if) 
and ip' are continuous at ±a; those Z that possess a decay eigenfunction we call decay 
eigenvalues or Gamow eigenvalues^ 

The remaining coefficients A±,B± are determined up to an overall factor by the 
requirement that both if) and its derivative if)' be continuous at ±a, the ends of the 
plateau. We will see that Imk < 0, so that, indeed, | "0(^)1 "~ > 00 exponentially as 
x — > ±oo. Continuity of if) and if)' determines the possible values of (not only A±,B± 
but also) Z and thus k, k. We have collected the details of the computations into 
Appendix [A] and report here the results. To describe the solutions, let 

Ap = - a = —V2rnAE = — . 41 

V2mAE A 

Ao is the de Broglie wavelength corresponding to the height AE of the potential plateau, 
and a is the width of the plateau in units of Aq. 



1 For Z = 0, the line for —a < x < a has to be replaced by Ao+Aix; for Z — —AE, if>(x) — D^+E-x 
for x < —a and if){x) = D + + E + x for x > a. 

2 Here is a look at the negative Z's that we excluded in this definition: In fact, for Z E (— oo, 0] \ 
{—AE} there exist no nonzero functions with C± = satisfying the eigenvalue equation ([55)1 for all 
x =/= ±a such that ip and ip' are continuous. (But we have not included the proof in this paper.) For 
Z = —AE, the coefficients C± are not defined, so the condition (j40|) makes no sense. 
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Proposition 1 If a > 10, then the number N of decay eigenvalues Z of fl33|) with 
\Z\ < AE/4 lies in the range a — 2 < N < a + 2. There is a natural way of numbering 
these eigenvalues as Z%, . . . , Z N . (There is no formula for the eigenvalue Z n , but it can 
be defined implicitly.) With each Z n is associated a unique (up to a factor) eigenfunction 
ip n , and \ip n (x)\ is exponentially increasing asi-> ±oo. Furthermore, forn <C a, 




Figure 8: Plot of |f/>n(:r)| 2 for an eigenfunction ip n with complex eigenvalue according 
to (|33|) with V(x) the plateau potential as in Figure [6j the parameters are n = 4 and 
a = 8; in units with a = 1, m = 1, and h~ = 1, this corresponds to A = 1/4 and 
AE = 32vr 2 = 315.8. 

The proof is given in Appendix |AJ the values ( 1421) are not (to our knowledge) in the 
literature so far. What can we read off about the lifetime r? In the regime a\/2mAE 3> 
h of large AE, which corresponds to a > 1, the n-th complex eigenvalue Z n with n<a 
is such that 

h 2 ir 2 n 2 2h 

ReZ n w , lmZ n « ReZ n . (43) 

8ma 2 aV2mAE 

(Readers familiar with the infinite well potential, V(x) = when —a < x < a and 
V(x) = oo when \x\ > a, which corresponds to the limit AE — > oo of very deep wells of 
the type shown in Figure 0, will notice that Re Z n given above actually coincides with 
the eigenvalues of the infinite well potential of length 2a.) Using (I3"6"l) and E = ReZ, 
one finds that 

V2mAE 1 [AE , , 

r = = — r cl = r qu , (44) 

the same value as specified in (I3"T1) . This completes our derivation of the lifetime fl3"TT) 
from complex eigenvalues. As we did for the potential step, we now also look at the 
question whether wave packets behave in the same way as the eigenfunctions, that is, 
whether a wave packet can remain on top of the plateau for the time span (13T|) . 

9 Wave Packets on the Plateau 

Now consider the long-time asymptotics for a normalized (square-integrable) wave func- 
tion of more or less definite energy such that the reflection coefficient is near 1, confined 



16 



initially in the plateau interval. As time proceeds, with the wave function slowly leaking 
out exponentially, a quasi-steady-state situation should be approached — in an expand- 
ing region surrounding the plateau — one that differs from a genuine steady state in that 
there is a global uniform overall exponential decay in time. (This state corresponds to 
a fixed point of the dynamics on projective space.) This situation should be described 
by a wave function of the sort discussed in the previous section, one that corresponds 
to a complex eigenvalue Z with negative imaginary part and that grows exponentially 
as \x\ — > oo. 

On the rigorous mathematical level, we can use the eigenf unctions with complex 
eigenvalue to establish that initial wave functions originating in the plateau region and 
involving only sufficiently low energy components have decay times given by r qu as in 
( 13"Tj) for the lowest relevant energy. As mentioned earlier, the complex Z n are not eigen- 
values of the quantum Hamiltonian in the standard mathematical terminology, and the 
corresponding eigenf unctions are not at all localized, but rather grow exponentially with 
increasing distance from the plateau. They nevertheless can be used for the mathemat- 
ical analysis. As a first step, we construct special states approximately localized in the 
plateau region which have the slow decay just described. This is the essence of Propo- 
sition [2] below. We then conclude in Corollary [T] that the slow decay applies to much 
more general initial states. 

Proposition 2 Let n be a fixed positive integer; keep a and m fixed and consider the 
regime h 2 /2ma 2 <C AE — > oo, so that a ^> n. There exists a normalized wave function 
4> n {x) such that (i) on the plateau <p n agrees (up to a factor) with the eigenfunction 
ip n with complex eigenvalue Z = Z n = E — ih/2r; (ii) the initial probability of the 
particle being on the plateau is close to 1; (Hi) for every < t < r, the time-evolved 
wave function <p n ^ = e~ lHt ^ h (j) n is close, on the plateau [—a, a], to (a factor times) the 
time-evolved eigenfunction ip n>t = e~ lZnt ^ n ij) n . Put differently, 

(i) 4> n (x) = A n ip n (x) for —a < x < a with some constant A n 6 C; 

(ii) l-/_ a j0 n (x)| 2 dx« 1; 

(™) I-a l^n,t(af) - A n ip n j(x)\ 2 dx < 1. 

The proof is included in Appendix [Bl As a consequence of Proposition [2], the amount 
of probability on the plateau indeed shrinks at rate 1/r, at least up to time r. In 
particular, the particle has probability ~ 1/e = 0.3679 to stay on the plateau until r. 

Let us describe how n is chosen. We make the n-th decay eigenfunction ip n normal- 
izable by multiplying it outside the plateau by a Gaussian function with width a > 0. 
We assume that a = 0(a), i.e., it does not grow with a. Finally, we multiply the re- 
sulting function by a constant A n which is so chosen that n has norm 1. Explicitly, 
then 




when x < —a 



when x > a, 

when — a < x < a. 



(45) 
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That is, (f) n is obtained from the eigenfunction ip n by "cutting off" the function around 
x = ±(a + a). Among the statements in Proposition [2J property (i) is immediate, and 
property (ii) is quite straightforward: Using the explicit form of A± and B± given in 
Appendix [A], one can compute that, in the regime n <^ a, 



\ip n (x)\ 2 dx ~ a . 



(46) 



while the integrals 



poo 




'a 





dx and 



x 



A 



dx 



(47) 



can be neglected in comparison, as they become small when a is large: They are of 
order 1/a 2 . This can be understood by noticing that |^ n (a)| 2 , the value at the left end 
of the interval [a, oo), is of order 1/a 2 , and the width a of the Gaussian is independent 



so 



of a, so that the integral is of order 1/a 2 , since \ip n \ 2 grows with rate of order 1/a 2 
slowly that the increase is negligible over an interval of length a. Property (iii) will be 
established in Appendix [0 



Corollary 1 For any initial wave function ip on the plateau with contributions only 
from eigenf unctions ip n with low n, i.e., 

"'max 

■ >^ c n ib n (x) when — a < x < a, 
rU): (48) 

otherwise, 

with a-independent n max and coefficients c n , the time-evolved wave function ip t = e~ lHt / n ip 
is close, on the plateau, to J2 n c n ip ntt , at least up to time min(ri, . . . , T„ max ), where 
r n = —h/2lmZ n . That is, 



/a 
\^t{x) - ^2c n ip ntt (x)\ 2 dx < 1 
-a 



(49) 



This means that any such wave packet ip will have a long decay time on the plateau, 
namely at least min(ri, . . . , T„ max ) (with each r n given by the quantum formula (T3TT) and 
not by the classical formula (jSDD Q; indeed, (Jl9"j) suggests that the decay time of ip is 
of the order of the largest r n with 1 < n < n max and significant |c n | 2 . The proof of 
Corollary [1] is simple; it is included in Appendix [Cl 

As another remark, the proof in Appendix [B] actually shows more than just Propo- 
sition [21 namely the following. Considering the eigenfunction ip n , we introduce the 
abbreviation 

v = — Refc, (50) 
m 
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which equals the speed at which an escaping particle moves away from the plateau. The 
maximal region in which we can expect a wave packet originally on the plateau to agree 
with an eigenfunction is the neighborhood of the plateau growing at speed v ; that is, the 
region at time t is the interval [—a — vt,a + vt]. And indeed, the wave function initially 
as in (I45p stays close to (a factor times) the eigenfunction in this growing region: 



for < t < r. Moreover, in the same way in which Corollary [T] is obtained from 
Proposition [2J one obtains the following from (1511) : For if) as in Corollary H) w = 
mm(v u . . . , w„ max ), v n given by ([50]) for if) n , and < t < min(r 1 , . . . , r nmax ), 



This justifies the qualitative description in the beginning of this section. 

As a final remark we note that the decay results described here, both qualitative and 
quantitative, presumably apply as well to the standard tunnelling situation in which a 
particle is confined inside a region by a potential barrier (a wall) that is high but not 
infinitely high, separating the inside from the outside. For this situation, more detailed 
results were obtained in [10J by other methods based on analytic continuation. 

10 Conclusions 

We have argued that paradoxical reflection and paradoxical confinement, the phe- 
nomenon that a quantum particle tends not to fall off a table, are real phenomena 
and not artifacts of the stationary analysis. We have pointed out that the effect is a 
robust prediction of the Schrodinger equation, as it persists when the potential step is 
not assumed to be rectangular but soft, and when the incoming wave is a packet of fi- 
nite width. We have provided numerical evidence and identified the relevant conditions 
on the parameters. We have explained why it is not a counter-argument to note that 
paradoxical reflection is classically impossible. We conclude that paradoxical reflection 
is a fact, not an artifact. Finally, we have computed that a state (of sufficiently low 
energy) on a potential plateau as in Figure [6] has a long decay time, no less than r qu 
given by (I3TI) . We conclude from this that a plateau potential can, for suitable param- 
eters, effectively be confining. Thus, the effect could indeed be used for constructing a 
(metastable) particle trap. 




a+vt 



(51) 




■a+vt 



(52) 
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A Solving the Plateau Eigenvalue Equation 

We now prove Proposition [TJ that is, we determine all decay eigenf unctions of fl33|) . as 
defined after (T4T)]) . The continuity of ip requires that 

A + e ika + A„e~ ika = B + e lka , (53) 
A + e- [ka + A_e ika = B_e lka , (54) 

and continuity of ip' that also 

k (A + e lka - A_e~ ika ) = kB + e ika , (55) 
k (A + e-' lka - A_e ika ) = -~kB_e ika . (56) 



Recall that both k and k can be complex. Since we assume AE > 0, and since, by (|33|) . 
k = k 2 + 2mAE/h 2 , we have that k + k ^ 0, and these equations are readily solved. 
First, we find the relations 

A _=^tlJl A (57) 
k + k 

Ik 

B ia(k-k) - A+ , (58) 

k + k 

k-k 

with the additional requirement that, since A + ^ for decay eigenf unctions, 

e i4afc . (60) 



2 

k + k 



k-k 

Let Ao and a be defined by (EH), and, in order to express k in natural units, let 

«:=^*. (61) 



Then 



and we have that 



k — k k — Vl + 
Thus fl60l) is equivalent to 



27T r 2lT i 

k = — — k , k = —Vl + K 2 , (62) 
Ao Aq 



k + h ' " Vl + '"" i • v'T~^) . (63) 



(k + v/TT^ 2 ") 4 = e i4 ™ a . (64) 
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The solutions of this equation coincide with those of the equation 



In 



(^k + \/T+ k 2 j = mna — i— 



(65) 



where n £ Z is arbitrary and In denotes the principal branch of the complex logarithm]^ 
Thus, with every decay eigenfunction if) is associated a solution k of fl65l) (with 
Reft > 0, since Re k > by definition of A;) and an integer n. Furthermore, n > — 1, 
because Re ft > and the imaginary part of the left hand side of f[6"5"j) must lie between 
— 7i and 7r. Conversely, with every solution k of (I6"5l) with ReK > there is associated 
a decay eigenvalue 

Z = k 2 AE (66) 

and an eigenfunction if) that is unique up to a factor: Indeed, (IBT!) and (|62|) provide the 
values of k and A; and imply (|66|) and (1601) ; Re/t > implies Z ^ (— oo,0], as well as 
Re A; > 0, so that indeed k = \ / 2mZ/h; k ±k ^ 0; A + can be chosen arbitrarily in 
C\ {0}, and if A_ and B± are chosen according to (I57|) - (I59|) then if) is nonzero (as, e.g., 
B + 7^ when k ^ and A + ^ 0) and a decay eigenfunction. We note that the condition 
ReK > is automatically satisfied when n > 2, as we can read off from the imaginary 
part of fl65l) using that In has imaginary part in (— tv, tt]. 

To determine ip explicitly, note that e i2afc |^ = (— l) n+1 , and thus = (-l) n+1 A + , 
B_ = (— l) n+1 B + ; setting A + = | and introducing the notation 



B := B + c 
we obtain that for odd n 



iak 



Aak 



K 



k + k 



\/l + K 2 + K 



(67) 



if)(x) = B x( x > a)e [k{x a) + \{ x < -a)e~ 



ik(x+a) 



and for even n 
ip(x) = B 



X (x > a)e iHx - a) - X (x < -a)e- fk{x+a) 



+ x(~ a ^ x < a) cos (kx) , (6J 



+ x(— a ^ x ^ a ) sm (k x ) ■ (69) 



We have used here the notation x(Q) to denote the characteristic function of a condition 
Q: 



x(Q) 



1 when Q is true, 
otherwise. 



(70) 



To sum up what we have so far, the decay eigenvalues are characterized, via (|66|) . 
through the solutions k of (1651) with Ren > 0. In order to study existence, uniqueness, 
and the asymptotics for a — > oo of these solutions, let us now assume, as in Proposition [fl 



3 For ( G C \ {0}, the equation e z = C has infinitely many solutions z, all of which have real part 
In and the imaginary parts of which differ by integer multiples of 2n; by ln^ we denote that z which 
has — 7r < Imz < 7r. 
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that a > 10 and \Z\ < AE/A. By virtue of fl66l) . the latter assumption is equivalent to 
\k\ < 1/2. We first show that solutions with \k\ < 1/2 must have \n\ < a + 2: Since In 
has imaginary part in (— 7r, it), flo^l) implies that Re k G ( ! |^, ! ^]) an d hence 

1 \n\ —2 

- > Id > |Re/d > — , (71) 

2 - i i - i 1 - 2a v ; 

or \n\ < a + 2. Next recall that for decay eigenvalues, n > — 1, so we obtain at this stage 
that the number of values that n can assume is at most a + 4, as the possible values are 
— 1, 0, 1, 2, . . . < a + 2. We will later exclude n = and n = —1. 

We now show that there exists a unique solution k of fl65l) for every n with |n| < a+2. 

Let 

= 7^- — In (k + VTT^) , (72) 

za 7TQ; V / 

so that fl65l) can equivalently be rewritten as the fixed point equation 

F(k) = k. (73) 

We use the Banach fixed point theorem [11] to conclude the existence and uniqueness 
of k. Since 

F'(k) = - — ^=L=, (74) 
Tea VI + K 

we have, by the triangle inequality, that 

= FTT 211/2 ^ 1 1 1 1211/2 • ( ?5 ) 

7ro;|l + k^I 1 /" 1 7TQ;|1 — l/tpl 1 /^ 

Let us consider for a moment, instead of \k\ < 1/2, the disk |k| < r for any radius 
< r < y/l - l/n 2 a 2 . There we have that \F'(k)\ < l/(naVl -r 2 ) =: K < 1. Thus, 
for any k, k' in the closed disk of radius r, \F(k') — F(n)\ < K\k' — k\, and, using 

\m\ = g, 

|F(k)| < |F(«) - F(0)| + |F(0)| <rf + ^<r, (76) 



2a 



provided that 



\n\ < 2ar{\ - K) . (77) 

Thus, in this case, F is a contraction in the ball of radius r, with a contraction constant 
of at most K. By the Banach fixed point theorem there is then a unique solution to the 
equation F(n) = k in the ball \k\ < r. Even though we are ultimately interested in the 
radius 1/2, let us set r = 1/a/2, which satisfies r < a/1 — l/7r 2 a 2 as a > 10; also ( 1771) 
is satisfied because |n | < a + 2 and a > 10 > 2(1 + 1/tt)/( v / 2 - 1) = 6.37. Hence, for 
every n with \n\ < a + 2, there is a unique solution n n with |/c n | < 1/V2. 

Getting back to the ball of radius 1/2, while some of the K n may have modulus greater 
than 1/2, we can at least conclude that there is at most one solution with modulus < 1/2 
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for every n with \n\ < a + 2. In addition, by setting r = 1/2, we obtain from f!77j) that 
\ K n\ < 1/2 for every n with |n| < a — 1. If n = 0, then -F(O) = and k = is the unique 
solution, which would lead to ip = 0. This excludes n = 0. Which of the solutions have 
Ren n > 0, as required for decay eigenvalues? For any n with \n\ < a + 2, let k„ be 
defined recursively by Kn = F{Kn^) with Kn = 0. Then, again by the Banach fixed 
point theorem for r = 1/ \/2, and 

For n = —1 and j = 1, this gives us that — < a~ 2 , and with k^ 1 ] = — l/2a 
and a > 10, we can conclude that Re/t_i < 0. This excludes n = —1. For n > 0, 
in contrast, the fact that \n n — Kn | < |n|a~ 2 allows us to conclude, with Kn = n/2a 
and a > 10, that ReK n > 0. Hence, the decay eigenvalues with \Z\ < AE/4 are in 
one-to-one correspondence with those K n , < n < a + 2, that have \k u \ < 1/2; the 
number of these K n must, as we have shown, be greater than a — 2 and less than or 
equal to a + 2. 

Furthermore, for these n n , lm.K n < 0: Computing K n explicitly yields 

4 2 ) = v - i— In f z/ + Vl + v 2 ) withz/ = — . (79) 
Tea \ J 2a 

Using (1781) as before, the claim follows if we can show that ImKn^ < —na~ 3 , which is 
equivalent to 

2tt 



— v < ln(u + Vl + v 2 ) . (80) 
a 

Note that < v < {a + 2) /2a < 0.6 since a > 10. In fact, (EDD holds for all < v < 1 
and a > 10 because \/l + v 2 > 1 and ln(z/ + l)/z/ for v > is a decreasing function^ 
whose value at u — 1 is In 2 = 0.693 > 0.628 = 27r/10 > 2-nja. As a consequence 
of Im/€ n < (and ReK n > 0), also Imk < 0, so that IVK^)! g rows exponentially as 
x — > ±oo. 

Now let us consider the asymptotics for n -C a. From (|78l we have that k„ is given 
by the right hand side of (1791) up to an error of order 0(na~ 3 ). Therefore, for integers 
n with < n < a we have that 

Tin . n ~ na . n 2 ^ 9 . _ n 2 AE ( ' . 2 \ . <N 

/c ~ i , fc^ i -, Z, n = k 2 AE w — l-i — . 81) 



2a 2aa a 4aa 2 ' " 4a 2 \ 7ra / 

B Derivation of the Lifetime Estimates for the Meta- 
stable States on the Plateau 

We now prove Proposition [2j What remains to be shown is property (iii); we will 
establish the stronger statement fl5"T|) . We will do this by constructing a function f(x, t) 



4 To see this, note that its derivative is of the form f{y)/g(v) with f{y) = v — (1 + v) ln(l + v) and 
g{v) = v 2 {l + v). Since g(v) > for v > 0, it suffices to show that f(v) < for v > 0. This follows 
from /(0) = and f'(y) = - hi(l + v) < for v > 0. 
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which remains close to the time-evolved eigenfunction in a growing region around the 
plateau. We then prove that this function forms an excellent approximation of (p n ^(x). 
f(x, t) will be defined by multiplying the time-evolved eigenfunction e~ itZn / h ip n again by 
Gaussians, but now with time-dependent parameters. 
To define f{x,t) we first introduce the abbreviation 



n 2 



Aaa 

and recall 



l3 = -lmk& — - (82) 



v = —Re k » — , (83) 
m ma 



whence k = j^v — i/3 with v, (3 > 0. We define further 

h 

R(t) = a + vt, b(t) = a 2 + i—t. (84) 

The Gaussians will be attached symmetrically to x — ±R(t) with a "variance" b(t), 
which yields explicitly 

( ±Be iH - x - a) -^ { - x ~ m) \ when x < -R(t), 

f{x,t) = A n e~ itz "/ h x J Be i~H*-a)-^(x-R(t))^ when x > R ^ ( 85 ) 

{^ n (x), when |x| < R(t). 

Clearly, f(x,0) = (j) n {x). Also f(x,t) is continuously differentiable in x for all t > 
because ip n is, and because the unnormalized Gaussian exp(— (x — /i) 2 /46) has, at its 
mean fx, value 1 and derivative 0. It is a short computation.] to check that for all t > 

(H - Z)f(x, t) = -£;dlf(x, t) + (V(x) - Z)f(x, t) 
= -&Mx- R(t),t) ± 9l (-x - R(t),t)] f{x, t), (86) 



2m 

\nJ<r — P(+ 

2m 

with 



»M = «(.>D(j|-^-^]. (87) 



In addition, we have 



ihd t f(x, t) = Zf(x, t)-£- [g 2 (x - R(t),t) ± g 2 (-x - R(t),t)} f(x, t), (88) 



5 The computation can be given the following mathematical justification: Since the potential V is 
bounded, by an application of the Kato-Rellich theorem [TH Theorem X.15], the Hamiltonian H = 
~'5m®x + ^ is self-adjoint on the domain of —9%. It can be easily checked that for any t the derivative 
d x f(x,i) is absolutely continuous in x, and thus the function f(-,t) belongs to the domain of H. This 
can be used to justify all the manipulations made here. Let us also use the opportunity to stress that, if 
we had not chosen the constants A± and B± in so that the function is continuously differentiable, 
then the addition of the Gaussian cut-off would have resulted in functions which are in L 2 (R) but which 
do not belong to the domain of H . Thus our estimates are not valid for such initial states. For more 
sophisticated mathematical methods to study such problems, see for instance |13) . 
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with g 2 = g± + g 3 where 



ga(y,t) = x(y > o) 



1 + 2/fy 
2b(t) ' 



(89) 



As for a fixed t, the function f(x,t) is square integrable, we can define a mapping 

t ^ F {t)^^f{.,t)-<t> n . (90) 

with F(t) G L 2 for all t > and F(0) =0. F is differentiable and by the above estimates 
for all t>0, 



d t F(t) = e 



itH/h 



h 



Hf(;t)+d t f(;t) 



jtH/h 



9(;t) 



where 



g(x, t) = i£- [g 3 (x - R(t),t) ± g 3 (-x - R(t),t)] f(x, t). 



(91) 



(92) 



As the derivative is continuous (in the L 2 -norm) in t, it can be integrated to yield 
F(t) = J* * ds d s F(s). Then, by the unitarity of the time evolution, we find 



||/(-,f)-M = \\F(t)\\ < fds \\d s F(s)\\ = fds \\g(; 

Jo Jo 



(93) 



Thus we only need to estimate the magnitude of J^ds \\g(-,s)\\. Let us first point 
out that for all \x\ > R(t), with y — \x\ — R(t), 



\f(x,t)\ 2 = |A l | 2 |S| 2 exp^2Re 



A z + ik{y + vt) -±bW y ' 

n 2 12 



(94) 



Here the argument of the exponential can be simplified using Z — ^—k — AE to 



a 



1 2 1 



A\b t \ 2y 2 Ct 2 



2£ 

Ct 



-y - ^ 



, with c t = 2/3\b(t)\a~ 1 . 



(95) 



Therefore, 



t] 1,2 = (i) 2 2 jo ^ lfiy + m ^ Uv ' t] |2 

2 

dy (1 + 2f3y) 2 exp(2Py - ^y 2 ) 
dx (1 + c 2 + c t x) 2 e~^ x2 
dx ((1 + c 2 ) 2 + c t V)e-^ 2 



2m|6*| / 2 



/i/3|A n ||B|Y Q io 

6 2 C t 



mact 



A(5 



~~ \ mac t J 4/3 



5 I\ C( 12 

ez * 



oo 



f gvgkygl ) I e Hv^((l + C 2 ) 2 + C 2 ) 

V 2ma / Cf 



(96) 
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For sufficiently large a and all < t < r m (2ma /Kim )a, 



c t <c T = 2(3\b(r)W- 1 = h\J^ + (^) 2 



r 2 



2 n 2 . ( a 2 \ 2 _ a 1 

V 4 + — a 2 < , (97) 



and therefore 



Since 



cr4aa 2 V \im 2 / na a 



| S MII<«. (« 



J(2(3a) 2 + (s(3hl(ma)) 2 > s— , (99) 

mcr 



we can estimate the integral over s by 



1 f , /mcr /mat 

ds ^ £ W^ 2 v^- < ioo > 



This proves that for all < t < t, 

\\f(.,t) - </>„,|| 2 < 4|A n | 2 |5| 2 — t « —-- « 1 . (101) 

mcr 7Tcr r a 

Since on the interval [—a — vt,a + vt], f(x, t) = A n ip n j(x), we have that 

t*a-\-vt roc 

\(j) njt {x) - A n ^ t (x)\ 2 dx < / \<f> nit (x) - f(x,t)\ 2 dx = \\f(;t) - (j) n , t \\ 2 < 1 , 
-a—vt J —oo 

(102) 

which is what we wanted to show. 

C Longevity of Wave Packets on the Plateau 

We now prove Corollary [TJ Define 

n=l 

and note that, since <p{x) = ip{x) for —a < x < a, and by virtue of property (ii) in 
Proposition [2j 

11-0 — 011 = ||0x(ic < —a or x > a) || < 5^ ~T~~ \\4>n x( x < — a or a; > a) || 1 . 
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As a consequence, also ip t will be close to 

0« = e- |K 



A r . 



so that, provided < t < r n for each n, 

(i>t ~ ^c n ^n,i)x(-a < x < a) I 

n 

(i>t - 0t)x(-« < x < a) + ^ - J^c^wW-a < ^ < a) 



< 1 



by virtue of property (iii) of Proposition [2J 



(103) 
(104) 
(105) 
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